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We employ Monte Carlo simulations to study a generalized three-dimensional complex \tlj\'^ theory 
of Ginzburg-Landau form and compare our numerical results with a recent quasi-analytical mean- 
field type approximation, which predicts first-order phase transitions in parts of the phase diagram. 
As we have shown earlier, this approximation does not apply to the standard formulation of the 
model. This motivated us to introduce a generalized Hamiltonian with an additional fugacity term 
controlling implicitly the vortex density. With this modification we find that the complex 
theory can, in fact, be tuned to undergo strong first-order phase transitions. The standard model 
is confirmed to exhibit continuous transitions which can be characterized by XY model exponents, 
as expected by universality arguments. A few remarks on the two-dimensional case are also made. 
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I. INTRODUCTION 

Since long the Ginzburg-Landau model has been con- 
sidered as paradigm for studying critical phenomena 
using field-theoretic techniques^ Perturbative calcula- 
tions of critical exponents and amplitude ratios of the 
Ising {n = 1), XY (n — 2), Heisenberg {n ~ 3) 
and other 0(n) spin models relied heavily on this field- 
theoretic formulation.^ Even though the spin models con- 
tain only directional fluctuations, while for ri-component 
Ginzburg-Landau fields with n > 2 directional and 
size fluctuations seem to be equally important, the two 
descriptions are completely equivalent, as is expected 
through the concept of universality and has been proved 
explicitly for superfluids with n — 2, where the spin 
model reduces to an XY modeL^ Therefore it appeared 
as a surprise when, on the basis of an approximate 
variational approach to the two-component Ginzburg- 
Landau model, Curty and Becki recently predicted for 
certain parameter ranges the possibility of first-order 
phase transitions induced by phase fluctuations. In sev- 
eral paperss^'§»i2iS this quasi-analytical^° prediction was 
tested by Monte Carlo simulations and, as the main re- 
sult, apparently confirmed numerically. If true, these 
findings would have an enormous impact on the theoreti- 
cal description of many related systems such as superfluid 
helium, superconductors, certain liquid crystals and pos- 
sibly even the electroweak standard model of elementary 
particle physicsi ^^d^ 

In view of these potential important implications for 
a broad variety of different fields we performed indepen- 
dent Monte Carlo simulations of the standard Ginzburg- 
Landau model in two and three dimensions in order to 
test whether the claim of phase-fluctuation induced first- 
order transitions is a real effect or not.^^ Our results 
clearly support the prevailing opinion that the nature 
of the transition is of second order. In turn this implies, 
of course, that the variational approximation employed 
in Ref. 4 is less reliable than originally thought in view of 
the apparent numerical confirmations. In onier_to_shed 
some light on the numerical results of Refs 



we 



generalized the standard model by adding a fugacity term 
which implicitly controls the vortex density of the model. 
The purpose of this paper is to present for this general- 
ized Ginzburg-Landau model results on its phase struc- 
ture as obtained from extensive Monte Carlo simulations. 
Employing finite-size scaling analyses we find numerical 
evidence that, by tuning the extra fugacity parameter, it 
is indeed possible to drive the system into a region with 
first-order phase transitions. 

The layout of the remainder of this paper is organized 
as follows. In Sec. we first recall the standard model, 
and then discuss its generalization and the observables 
used to map out the phase diagram. Next we describe the 
employed simulation techniques in Sec. IIIII The results 
of our simulations are presented in Sec. IIVI where we 
first discuss the three-dimensional case in some detail and 
then add a few brief comments on the two-dimensional 
model to complete the physical picture. Finally, in Sec.lVl 
we conclude with a summary of our main findings. 



II. MODEL AND OBSERVABLES 

The standard complex or two-component Ginzburg- 
Landau theory is defined by the Hamiltonian 



7>0, (1) 



where ^(r) = + i^yi^) = \^{r)\e^'^^'^^ is a com- 

plex field, and a, b and 7 are temperature independent 
coefficients derived from a microscopic model. In order 
to carry out Monte Carlo simulations we put the model 
on a d-dimensional hypercubic lattice with spacing a. 
Adopting the n otation of Ref.0, we introduce scaled vari- 
ables ^ = tp/ y/ {\a\/b) and u = r/^, where ^ — vt/H 
is the mean-field correlation length at zero temperature. 
This leads to the normalized lattice Hamiltonian 
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where /i denotes the unit vectors along the d coordinate 
axes, N = L'^ IS the total number of sites, and an unim- 
portant constant term has been removed. The param- 
eter Vo merely sets the temperature scale and can thus 
be absorbed in the definition of the reduced temperature 
f = T/Vo. 

After these rescalings, and omitting the tilde on -0, cr, 
and T for notational simplicity in the rest of the paper, 
the partition function Z considered in the simulations is 
then given by 
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and / Dip Dip = J D Keijj D Imip stands short for inte- 
grating over all possible complex field configurations. 

In Ref. Il.'j we have shown, that the disagreement men- 
tioned above is caused by an incorrect sampling of the 
Jacobian which emerges from the complex measure in 
when transforming the field representation to polar 
coordinates, i/'n — -Rn(cos((/)„), sin((/)„)). When updating 
in the simulations the modulus i?„ — Itpnl and the an- 
gle (/)„, one has to rewrite the measure of the partition 
function as 



Z = 



where DR = Yln=i '^^n and R = Y\n=i Rn is the Jaco- 
bian of this transformation. While mathematically in- 
deed trivial (and of course properly taken into account 
m Ref.!?), this fact may easily be overlooked when cod- 
ing the update proposals for the modulus and angle in 
a Monte Carlo simulation program. While for the an- 
gles it is correct to use update proposals of the form 
0n ~* <t>n + 54) with — A0 < <50 < A0 (where A</) is 
chosen such as to assure an optimal acceptance ratio), 
the similar procedure for the modulus, i?„ Rn + 5R 
with —AR < SR < AR, would be incorrect since this ig- 
nores the Rn factor coming from the Jacobian. In fact, if 
we purposely ignore the Jacobian and simulate the model 
© (erroneously) without the i?- factor, then we obtain a 
completely different behavior than in the correct case, cf. 
e.g. Fig. 12 below. As already mentioned above these re- 
sults reproducei^ those in Refs. |^ and 0, and from this 
data one would indeed conclude evidence for a first-order 
phase transition when a is small. With the correct mea- 
sure, on the other hand, we have checked that no first- 
order signal shows up down to a = 0.01. 
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To treat the measure in Eq. JBJ properly one can ei- 
ther use the identity RndRn — dR^/2 and update the 
squared moduli i?^ = IV'nP according to a uniform mea- 
sure (where the update proposal Rn R^ + S with 
— A < (5 < A is correct), or one can introduce an ef- 
fective Hamiltonian, 
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with K = 1 and work directly with a uniform measure 
for Rn- The incorrect omission of the i?-factor in © 
is equivalent to setting k ~ 0. It is well knowrkii that 
the nodes Rn = correspond to core regions of vor- 
tices in the dual formulation of the model. The Jacobian 
factor R (or equivalently the term — ^lni?„ in -ffoff) 
tends to suppress field configurations with many nodes 
Rn = 0. If the i?-factor is omitted, the number of nodes 
and hence vortices is relatively enhanced. It is thus at 
least qualitatively plausible that in this case a discon- 
tinuous, first-order "freezing transition" from a vortex 
dominated phase can occur, as is suggested by a similar 
mechanism for the XY modelii4iSii& and defect-models 
of meltingiii^. 

In the limit of a large parameter a, it is easy to read 
off from Eq. that the modulus of the field is squeezed 
onto unity and once hence expects that irrespectively of 
the value of k the XY model limit is approached with its 
well-known continuous phase transition in three dimen- 
sions (3D) at Tc w 2.2 respectively Kosterlitz-Thouless 
(KT) transition in two dimensions (2D) at Tkt ~ 0.9. 
While for the standard model with k = 1, this behavior 
should qualitatively persist for all values of cr, from the 
numerical results discussed above one expects that for 
K = the order of the transition turns first-order below 
a certain (tricritical) cr-value. The purpose of this pa- 
per is to elucidate this behavior further by studying the 
phase diagram in the tr-K-plane, i.e., by considering an 
interpolating model with k varying continuously between 
and 1. 

To be precise we always worked with the proper func- 
tional measure in Eq. © and replaced the standard 
Hamiltonian H by 
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where we have introduced the parameter S = 1 — k, 
such that S = {k = 1) corresponds to the standard 
model and i5 = 1 (k = 0) to the previously studied mod- 
ified model with its first-order phase transition for small 
enough a. 

In order to map out the phase diagram in the a-K- 
respectively cr-i5-plane, we have measured in our simula- 
tions to be described in detail in the next section among 
other quantities the energy density e = {H)/N, the spe- 
cific heat per site Ci, = {{H^) — {H)^)/N, and in particular 
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the mean-square amplitude 
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which will serve as the most relevant quantity for compar- 
ison with previous workiS*SiSia. For further comparison 
and in order to determine the critical temperature, the 
helicity modulus, 
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was also computed. Notice that the helicity modulus 
is a direct measure of the phase correlations in the 
direction of /i. Because of cubic symmetry all direc- 
tions /i are equivalent, and we always quote the average 
r = X])l=i ^M- infinite-volume limit, F is 

zero above and different from zero below Tc. We also 
have measured the vortex density v (of vortex points in 
2D and vortex lines in 3D). The standard procedure to 
calculate the vorticity on each plaquette is by considering 
the quantity 
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where . . . , 4)4 are the phases at the corners of a pla- 
quette labeled, say, according to the right-hand rule, 
and [a]2-K stands for a modulo 27r: [a]2TT — a + 27rrt, 
with n an integer such that a + 2Tm G (— 7r,7r], hence 
771 = 77 1 2 + + 7734 + 7*41 . If 7)7, ^ 0, there exists a vortex 
which is assigned to the object dual to the given plaque- 
tte (a site in 2D and a link in 3D). Hence, in two dimen- 
sions, *777, the dual of m, is assigned to the center of the 
original plaquette. In three dimensions, the topological 
point charges are replaced by (oriented) line elements *li 
which combine to form closed networks ( "vortex loops" ) . 
The vortex "charges" *m or *li can take three values: 
0, ±1 (the values ±2 have a negligible probability). The 
quantities 
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serve as a measure of the vortex density. We further 
analyzed the Binder cumulant, 
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where fl = (/ia;,/ij,) with 
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FIG. 1: Mean-square amplitude of the standard three- 
dimensional complex Ginzburg-Landau model with n — 1 and 
a = 0.25 on a 10'' cubic lattice. 



is the magnetization per lattice site of a given configura- 
tion. 



III. SIMULATION TECHNIQUES 

Let us now turn to the description of the Monte Carlo 
update procedures used by us. To be on safe grounds, 
we started with the most straightforward (but most in- 
efficient) algorithm known since the early days of Monte 
Carlo simulations: The standard Metropolis algorithm^Si. 
Here the complex field tpn is decomposed into its Carte- 
sian components, 4'n = ipx,n + ii^y,n- For each lattice 
site a random update proposal for the two components 
is made, e.g. Va;,™ ^ '^x,n + Sipx,n with JVx.n S [-A, A], 
and in the standard fashion accepted or rejected accord- 
ing to the energy change SHg^n- The parameter A is usu- 
ally chosen such as to give an acceptance rate of about 
50%, but other choices are permissible and may even re- 
sult in a better performance of the algorithm (in terms 
of autocorrelation times). All this is standard'^^ and 
guarantees in a straightforward manner that the complex 
measure DipDip in the partition function Q is treated 
properly. 

The well-known drawback of this algorithm is its crit- 
ical slowing down (large autocorrelation times) in the 
vicinity of a continuous phase transition^S, leading to 
large statistical errors for a fixed computer budget. To 
improve the accuracy of our data we therefore employed 
the single-cluster algorithm.'^^ to update the direction of 
the field^, similar to simulations of the XY spin model^— . 
The modulus of ip is updated again with a Metropolis 
algorithm. Here some care is necessary to treat the mea- 
sure in ^ properly (see above comments). Per mea- 
surement we performed one sweep with the Metropolis 
algorithm and 77 single-cluster updates. For all simula- 
tions in two and three dimensions the number of cluster 
updates was chosen such that 77(|C|) ~ L'^ = N, where 
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(|C|) is the average cluster size. Since (|C|) scales with 
system size as the susceptibility, x = ^{fl^) — L^^" ^ and 
7/j/ = 2 — 77 = 7/4at the Kosterlitz-Thouless transition 
in 2D and 7/j/ = 2 — 77R:!2in 3D, n was chosen oc L^^^ 
in 2D and oc L in 3D. In the 2D case most of the sim- 
ulations were performed for L = 10, 20, and 40, and in 
3D we usually studied the lattice sizes L = 10, 15, 20, 
and 30. For each simulation point we thermalized with 
500 to 1 ODD sweeps and averaged the measurements over 
10 000 sweeps. In the cases of strong first-order phase 
transitions we employed a variant of the multicanonical 
scheme^* where the histogram of the mean modulus is 
flattened instead that of the energy. All error bars are 
computed with the Jackknife method^^. In the following 
we only show the more extensive and accurate data set 
of the cluster simulations, but we tested in many repre- 
sentative cases that the Metropolis simulations coincide 
within error bars, for an example see Fig.^ 



IV. RESULTS 

A. Three dimensions 

In the first set of simulations we concentrated on the 
two most characteristic cases k = and k = 1 and per- 
formed temperature scans on a 15^ lattice for various val- 
ues of the parameter a. Our results for the mean-square 
amplitude, the helicity modulus and the vortex-line den- 
sity are compared for the two cases in Fig. |21 In the 
plots for K = on the left side, we see that all three 
quantities exhibit quite pronounced jumps for small a- 
values, which is a clear indication that in this regime the 
phase transition is of first order. At a = 0.25, for ex- 
ample, we observe already on very small lattices a clear 
double-peak structure for the distributions of the energy 
and mean-square amplitude as well as the mean modulus 
IV'I = jj^n=i I'/'" I ■^^hich is depicted in Fig. El Notice 
that already for the extremely small lattice size of 4"^ 
the minimum between the two peaks is suppressed by 
more than 20 orders of magnitude. This is an unam- 
biguous indication for two coexisting phases and thus 
clearly implies that the model undergoes a first-order 
phase transition in the small cr- regime for n = 0. Due to 
the pronounced metastability these simulations had to be 
performed with a variant of the multicanonical scheme^* 
where, instead of flattening the energy histogram, extra 
weight factors for the mean modulus were introduced. 
With this simulation technique we overcome the diffi- 
culty of sampling the extremely rare events between the 
two peaks of the canonical distribution. A closer look at 
the K = plots shows that the crossover from second- 
to first-order transitions happens around at ~ 2.5. For 
the standard model with k = 1, on the other hand, we 
observe for all a-values a smooth behavior, suggesting 
that the XY model like continuous transition persists 
also for small cr-values. This is clearly supported by a 
single-peak structure of all distributions just mentioned. 



for the case of the mean modulus see Fig. O This sup- 
ports the prevailing opinion that the standard complex 
■01'' model always undergoes a second-order phase tran- 
sition. In fact, we have checked that down to cr = 0.01 no 
signal of a first-order transition can be detected for the 
standard model parameterized by k = 1. The resulting 
transition lines in the cr- T-plane for k = and k = 1 are 
sketched in Fig. 01 with the thick line for k = indicating 
the approximate regime of first-order phase transitions. 

Next we concentrated on the small cr regime and per- 
formed a rough finite-size scaling (FSS) analysis for cr = 
0.25 on moderately large 10^, 15^, 20^, and 30^ lattices. 
In Fig.[5lwe compare results for the energy, mean-square 
amplitude © , helicity modulus H1Q(I and vortex-line den- 
sity (|13|l for K = and k = 1. Apart from the transi- 
tion region where a strong size dependence is of course 
expected, we notice only a small dependence on the vari- 
ation of the lattice size. On the basis of these results, 
we do not expect a significant change of the qualitative 
behavior for much larger lattices and hence used similar 
moderate lattice sizes for most of our further investiga- 
tions. 

To exemplify the big differences between the models 
with K = and k = 1, we choose in the following the 
case cr — 1.5, where we shall characterize for both re- 
values the phase transitions in some detail. Let us start 
with the non-standard case k — 0, where the first-order 
phase transition around T « 0.36 is also pronounced but 
much less strong than for cr = 0.25. Still, in order to get 
sufficiently accurate equilibrium results, the simulations 
for lattices of size L = 4,6,8,10,12,14,15, and 16 had 
to be performed again with our modulus variant of the 
multicanonical method. As can be inspected in the his- 
togram plots for the mean modulus shown in Fig. El the 
frequency of the rare events between the two peaks in the 
canonical ensemble for a 16^ lattice is about 50 orders of 
magnitude smaller than for configurations contributing 
to the two peaks. 

In order to characterize the transition more quantita- 
tively we estimated the interface tension^^^ 
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where p™ax 



is the value of the two peaks and P™'" de- 
notes the minimum in between. Here we have assumed 
that for each lattice size the temperature was chosen 
such that the two peaks are of equal height which can 
be achieved by histogram reweighting. The thus de- 
fined temperatures approach the infinite-volume transi- 
tion temperature as 1/i'*, and for the final estimate of 
= limi^oo F[, we performed a fit according to^ 



F? = F" 
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As is shown in Fig. |B1 the finite-lattice estimates Fl are 
clearly nonzero. The infinite-volume extrapolation H17|l 
tends to increase with system size and yields a compara- 
bly large interface tension of F" = 0.271(5). 
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FIG. 2: Mean-square amplitude (IV'P)) helicity modulus and vortex-line density {v) of the three-dimensional generalized 
complex Ginzburg-Landau model on a 15"^ cubic lattice for different values of the parameter a = 0.25, . . . , 3.0 for the case k = 
(left) and the standard formulation with k = 1 (right). 



Let us now turn to the second generic case, k — 
1, where the model definitely exhibits for a ^ 1.5 
a second-order phase transition around f3 = 1/T « 
0.8. To confirm the expected critical exponents of 
the 0(2) or XY model universality class, we simulated 
here close to criticality somewhat larger lattices of size 
L = 4,8,12,16,20,24,32,40, and 48 and performed a 
standard FSS analysis. From short runs we first es- 



timated the location of the phase transition to be at 
/3o = 0.7795 « (3c- In the long runs at (3o we recorded 
the time series of the energy density e = E/N, the 
magnetization /2, the mean modulus and the mean- 
square amplitudes^ as well as the helicity mod- 
ulus and the vorticity v. After an initial equili- 
bration time we took about 1 000 000 measurements for 
each lattice size. Applying the reweighting technique 
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FIG. 3: Top: Histogram of the mean modulus \tp\ on a log- 
arithmic scale for a 4^ cubic lattice, k = and a = 0.25, 
reweighted to the temperature To ~ 0.0572 where the two 
peaks are of equal height. Bottom: Histogram for the same 
quantity and lattice size at T — 1.1 close to the second-order 
phase transition for n = 1 and a = 0.25. 

we first determined the maxima of the susceptibihty, 
X' = A^((/2^)-(|/i|)^), o{d{\fl\)/d/3, and of the logarithmic 
derivatives dln{\fl\) /df3 and d\\i{fp) / dfi . The locations of 
these maxima provide us with four sequences of pseudo- 
transition points /3max(i) for which the scaling variable 
X = {/3max(i) — Pc)L^'^ should be constant. Using this 
fact we then have several possibilities to extract the crit- 
ical exponent v from (linear) least-squares fits of the FSS 
ansatz dUL/d(3 =^ L^'" fo{x) or dln(|/2|P)/d/3 ^ L^'" fp{x) 
to the data at the various Anax(-^) sequences. The qual- 
ity of our data and the fits starting at Lmin — 8, with 
goodness-of-fit parameters Q — 0.85 — 0.90, can be in- 
spected in Fig. [7| All resulting exponent estimates and 
consequently also their weighted average, 

= 1.493(7), J/ = 0.670(3), (18) 

are in perfect agreement with recent high-precision 
Monte Carlo estimates for the XY model universality 
classjSSiSS Note that hyperscaling implies a = 2 — 2>v = 
—0.010(9), which also favorably compares with recent 




0.5 1 1.5 2 2.5 3 





FIG. 4: Transition lines in the cr-T-plane for k = and k = 1. 
The thick line for k = indicates first-order phase transitions 
while all other transitions are continuous. 

spacelab experiments on the lambda transition in liquid 
helium iffi 

Assuming thus X/v — 1.493 we can improve our esti- 
mate for from linear least-squares fits to the scaling 
behavior of the various /3max sequences. The combined 
estimate from the four sequences is /3c = 0.780 08(4). To 
extract the critical exponent ratio ^/v we can now use 
the scaling relation for the susceptibility x — A^(m^) — 
a + bL''/'' at Pc- For L > 16 we obtain from a FSS fit 
with Q = 0.70 the estimate of 

= 1.962(12)[9], (19) 

where we also take into account the uncertainty in our es- 
timate of (3c] this error is estimated by repeating the fit at 
/9c± A/3c and indicated by the number in square brackets. 
Here we find a slight dependence of this value on the lower 
bound of the fit range [I/min,48], i.e., one would have 
to include larger lattices for a high-precision estimate of 
the critical exponent ratio 7/1^, but this was not our ob- 
jective here. Still, these results are in good agreement 
with recent high-precision estimates in the literatureSSiSS 
and clearly confirm the expected second-order nature of 
the phase transition in the standard complex \ip\'^ model, 
governed by XY model critical exponents. 

A similar set of simulations at a = 0.25 for lattice 
sizes L = 4, 8, 12, 14, 16, 20, 24, 28, 32, and 40 gave the 
exponent estimates l/v — 1.498(9), v = 0.668(4) and 
= 1.918(71) [8] (at (i^ = 0.9284(4)), which are less 
accurate but again compatible with the XY model uni- 
versality class. At any rate these results definitely rule 
out the possibility of a first-order phase transition in 
the standard model at small cr-values. When going to 
even smaller cr-values, the FSS analysis is more and 
more severely hampered by the vicinity of the Gaussian 
fixed point which induces strong crossover scaling effects. 
Since consequently very large system sizes would be re- 
quired to see the true, asymptotic (XY model like) crit- 
ical behavior we have not further pursued our attempts 
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FIG. 5: Energy density e, mean-square amplitude {|V'P)i helicity modulus and vortex-line density v on 10^,15^,20^ and 
30^ cubic lattices for a = 0.25 and k = (left) respectively k = 1 (right). 
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FIG. 6: Top: Histogram of the mean modulus \^\ for k = 
and o" = 1.5 on a logarithmic scale for various lattice sizes 
ranging from L = 4 (top curve) to L = 16 (bottom curve), 
reweighted to temperatures where the two peaks are of equal 
height. Bottom: FSS extrapolation for L > 6 of the interface 
tension _F£, yielding the infinite- volume limit F" = 0.271(5). 



in this direction. Here we only add the remark that for 
a — 0.01 the energy and magnetization distributions ex- 
hibit a clear single-peak structure for all considered lat- 
tice sizes up to L = 20, showing that in the standard 
model with k = 1 a phase-fluctuation induced first-order 
phase transition is very unlikely even for very small a 
values. 

We also checked the critical behavior along the line of 
second-order transitions for k — 0. Specifically, at tr = 5, 
i.e., sufficiently far away from the crossover to first-order 
transitions at crt « 2.5, we obtained from FSS fits to data 
for lattices of size L = 4, 8, 12, 16, 20, 24, 28, 32, and 40 
the exponent estimates l/u — 1.489(7), — 0.671(3) and 
= 1.913(82) [13] (at Pc = 0.97253(4)). As expected 
by symmetry arguments, also these results for the second- 
order regime of the k = variant of the model are in 
accord with the XY model universality class. 

In a second set of simulations we explored the two- 
dimensional a-K parameter space of the generalized 
Ginzburg-Landau model in the orthogonal direction by 
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FIG. 7: Least-square fits for k — 1 and cr = 1.5 on a log- log 
scale, using the FSS ansatz df{iJ,)/dj3 oc L^'''^ at the maxima 
locations. The fits using the data for L > 8 lead to an overall 
critical exponent 1/u = 1.493(7) or u = 0.670(3). 
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FIG. 8: Log-log plot of the FSS of the susceptibility for k = 1 
and a = 1.5 at ^ = 0.780 08 « /3c. The line shows the three- 
parameter fit a -I- hU'^" , yielding for L > 16 the estimate 
-l/u = 1.962(12). 



performing simulations at fixed cr values and k varying 
from K = to 1. For most cr- values we concentrated on 
the crossover region between first- and second-order tran- 
sitions when varying n. For two selected values, a — 0.25 
and cr = 1.5, we studied the k dependence more system- 
atically by simulating all values from k = to 1 in steps 
of 0.1. In addition we performed two further runs in the 
crossover regime at k = 0.85 and 0.95 for cr = 0.25 as well 
as at K = 0.15 and k, = 0.25 for cr — 1.5. In Fig.|51we show 
the resulting mean-square amplitudes for all simulated 
values of K at CT = 0.25 as a function of the temperature, 
indicating again that for small k, the transitions are first- 
order like while for k, closer to unity the expected second- 
order transitions emerge. From Fig. |51 we read off that 
for cr = 0.25 the crossover between the two types of phase 
transitions happens around Kt(cr = 0.25) « 0.8, and the 
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FIG. 9: The k dependence of the mean-square amphtude 
(l-i/jp) as a function of temperature on a 15'^ lattice for a = 
0.25. 
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FIG. 10: Phase diagram in the K-T-plane of the three- 
dimensional generalized complex Ginzburg-Landau model for 
a = 0.25 and a = 1.5. The transitions along the thick line for 
K < Kt are of first order, and the transitions for k > Kt are 
of second order. The points labeled hit at the intersection of 
these two regimes are tricritical points. 



analogous analysis for a = 1.5 yields Kt{<J — 1.5) « 0.2. 
The resulting transition lines for these two cr-values are 
plotted in Fig. ^| where the thick lines indicate again 
first-order phase transitions. 

Finally, by combining all numerical evidences collected 
so far with additional data not described here in detail, 
we find the phase structure in the cr-K-plane depicted in 
Fig. 111! All points in the lower left corner for small a and 
small K exhibit temperature driven first-order phase tran- 
sition when the temperature is varied, while all points in 
the upper right corner display a continuous transition of 
the XY model type. This means in particular that for 
the standard model parameterized by k = 1 this is al- 
ways true. Quantitatively the XY model is reached for 
all K-values in the limiting case a — > oo. 



FIG. 11: Phase structure in the a-K-plane of the generalized 
complex Ginzburg-Landau model in three dimensions, sepa- 
rating regions with first- and second-order phase transitions, 
respectively, when the temperature is varied. All continuous 
transitions fall into the universality class of the XY model 
which is approached for all k- values in the limit a — > oo. 

B. Two dimensions 

We conclude the paper with a few very brief remarks 
on the two-dimensional generalized model where the 
Kosterlitz-Thouless nature of the standard XY model 
transition would require more care for a precise study. 
Here we only report results of some runs at ct = 1 for 
10^, 20^, and 40^ square lattices. As the main result, 
we find that the standard observables e, (IV'l^)i T, and v 
exhibit qualitatively the same pattern as in three dimen- 
sions. This is demonstrated in Fig. E| where again the 
two cases k = and k = 1 are compared. For k = 0, 
the data are indicative of a first-order transition around 
T K, 0.2, while the behavior of the standard model with 
K = 1 is consistent with the expected Kosterlitz-Thouless 
transition around T k, 0.4. Note in particular that (only) 
the data for k — 1 are compatible with the expected uni- 
versal jump of the helicity modulus at Tc, Ty = (2/7r)T, 
indicated by the straight line in the corresponding plots. 
A careful investigation of the first-order transitions in the 
generalized model with k = will be reported elsewhere. 

V. SUMMARY 

The possibility of a phase-fluctuation induced first- 
order phase transition in the standard three-dimensional 
Ginzburg-Landau model as suggested by approximate 
variational calculations^ cannot be confirmed by our nu- 
merical simulations down to very small values of the pa- 
rameter a. Our results suggest, however, that a gener- 
alized Ginzburg-Landau model can be tuned to undergo 
first-order transitions by a mechanism similar to that dis- 
cussed in Ref. when varying the parameter k of an 
additional ^ In i?„ term in the generalized Hamiltonian 
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FIG. 12: Energy density e, mean-square amplitude (IV'P), helicity modulus and vortex density v of the two-dimensional 
model on 10^20^ and 40^ square lattices for a — 1 and k = (left) respectively k = 1 (right). The straight lino in the F^ 
plots indicate the universal KT jump F^^ = (2/7r)T at T = Tc, which clearly is only compatible with the data for the standard 
model with k = 1. 



11 



As in Ref. this can be understood by a duality 
argument. For < k < 1 the extra term reduces the ra- 
tio of core energies of vortex Hues of vorticity two versus 
those of vorticity one, and this leads to the same type of 
transition as observed in defect melting of crystals. 

The phase transitions of the standard model as well as 
the continuous transitions of the generalized model are 
confirmed to be governed by the critical exponents of the 
XY model or 0(2) universality class, as expected by gen- 
eral symmetry arguments. For the generalized model it 
would be interesting to analyze in more detail the tricrit- 
ical points separating the regions with first- and second- 
order phase transitions. Such a study, however, is quite 
a challenging project and hence left for the future. 

Exploratory simulations of the two-dimensional case. 



where the standard model exhibits Kosterlitz-Thouless 
transitions, indicate that a similar mechanism can drive 
the transition of the generalized model to first order also 
there. 
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